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Abstract — The goal of this paper is to establish a general 
approach for analyzing queueing models with repeated inho- 
mogeneous vacations. The server goes on for a vacation if the 
inactivity prolongs more than the vacation trigger duration. 
Once the system enters in vacation mode, it may continue for 
several consecutive vacations. At the end of a vacation, the server 
goes on another vacation, possibly with a different probability 
distribution; if during the previous vacation there have been 
no arrivals. However the system enters in vacation mode only 
if the inactivity is persisted beyond denned trigger duration. 
I In order to get an insight on the influence of parameters on 
the performance, we choose to study a simple M/G/l queue 
(Poisson arrivals and general independent service times) which 
has the advantage of being tractable analytically. The theoretical 
model is applied to the problem of power saving for mobile 
devices in which the sleep durations of a device correspond to 
the vacations of the server. Various system performance metrics 
such as the frame response time and the economy of energy are 
derived. A constrained optimization problem is formulated to 
maximize the economy of energy achieved in power save mode, 
with constraints as QoS conditions to be met. An illustration of 
the proposed methods is shown with a WiMAX system scenario 
to obtain design parameters for better performance. Our analysis 
allows us not only to optimize the system parameters for a given 
traffic intensity but also to propose parameters that provide the 
best performance under worst case conditions. 

Index Terms — M/G/l queue with repeated vacations, power 
save mode, system response time, gain optimization. 

I. Introduction 

Power save/sleep mode operation is the key point for energy 
efficient usage of mobile devices driven by limited battery 
lifetime. Current standards of Mobile communication such as 
WiFi, 3G and WiMAX have provisions to operate the mobile 
station in power save mode in case of low uses scenarios. A 
mobile operating in power save or sleep mode saves the battery 
energy and enhances lifetime but it also introduces unwanted 
delay in serving data packets arriving during sleep duration. 
Though energy is a major aspect for handheld devices, delays 
may also be crucial for various QoS services such as voice 
and video traffic. Mobility extension of WiMAX JT] is one 
of the most recent technologies whose sleep mode operation 
is being discussed in detail and standardized. 

The IEEE 802. 16e standard [JTJ defines 3 types of power 
saving classes. 

• Type I classes are recommended for connections of Best- 
Effort (BE) and Non-Real Time Variable Rate (NRT-VR) 
traffic. Under the sleep mode operation, sleep and listen 
windows are interleaved as long as there is no downlink 
traffic destined to the node. During listen windows, the 
node checks with the base station whether there is any 
buffered downlink traffic destined to it in which case 

e-mail: amar. azad@ieee.org 



it leaves the sleep mode. Each sleep window is twice 
the size of the previous one but it is not greater than 
a specified final value. A node may awaken in a sleep 
window if it has uplink traffic to transmit. 

« Type II classes are recommended for connections of 
Unsolicited Grant Service (UGS) and Real-Time Variable 
Rate (RT-VR) traffic. All sleep windows are of the same 
size as the initial window. Sleep and listen windows are 
interleaved as in type I classes. However, unlike type I 
classes, a node may send or receive traffic during listen 
windows if the requests handling time is short enough. 

• Type III classes are recommended for multicast con- 
nections and management operations. There is only one 
sleep window whose size is the specified final value. 
At the expiration of this window, the node awakens 
automatically. 

The related operational parameters including the initial and 
maximum sleep window sizes can be negotiated between the 
mobile node and the base station. 

The sleep mode operation of IEEE 802. 16e, more specifi- 
cally the type I power saving class, has received an increased 
attention recently. In 0, the base station queue is seen as 
an M/GI/l/N queueing system with multiple vacations; an 
embedded Markov chain models the successive (increasing in 
size) sleep windows. Solving for the stationary distribution, 
the dropping probability and the mean waiting time of down- 
link packets are computed. Analytical models for evaluating 
the performance in terms of energy consumption and frame 
response time are proposed in Q, H and supported by 
simulation results. While @ considers incoming traffic solely, 
both incoming and outgoing traffic are considered in 0. 
In 13, the authors evaluate the performance of the type I power 
saving class of IEEE 802. 16e in terms of packet delay and 
power consumption through the analysis of a semi-Markov 
chain. 

Power save mode in systems other than the IEEE 802. 16e 
have also been studied; hereafter we cite some of these studies. 
In 16], the authors evaluate the energy consumption of various 
access protocols for wireless infrastructure networks. The 
sleep mode operation of Cellular Digital Packet Data (CDPD) 
has been investigated through simulations in Q and analyti- 
cally in HI. To efficiently support short-lived sessions such as 
web traffic, a bounded slowdown method - that is similar to 
type I power saving classes in the IEEE 802. 16e - is proposed 
for the IEEE 802.11 protocol in (9j. Last, the power saving 
mechanism for the 3G Universal Mobile Telecommunications 
System (UMTS) system is evaluated in iflOl . 

In this paper, we propose a queueing-based modeling frame- 
work that is general enough to study many of the power 
save operations described in standards and in the literature. 
In particular, our model enables the characterization of the 
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performance of type I and type II power saving classes as de- 
fined in the IEEE 802. 16e standard UJ. The system composed 
of the base station, the wireless channel and the mobile node is 
modeled as an M/G/l queue with repeated inhomogeneous 
vacations. Traffic destined to the mobile node awaits in the 
base station as long as the node is in power save mode. 
When the node awakens, the awaiting requests start being 
served on a first-come-first-served basis. The service consists 
of the handling of a frame at the base station, its successful 
transmission over the wireless channel and its handling at 
the node. Analytical expressions for the distribution and/or 
the expectation of many performance metrics are derived 
yielding the expected frame transfer time and the expected 
gain in energy. We formulate an optimization problem so 
as to maximize the energy efficiency gain, constrained to 
meeting some QoS requirements. We illustrate the proposed 
optimization scheme through four application scenarios. 

Although we have motivated our modeling framework using 
power saving operation in wireless technologies, it is useful 
whenever the system can be modeled by a server with repeated 
vacations. The structure of the idle period is general enough 
to accommodate a large variety of scenarios. 

There has been a rich literature on queues with vacations, 
see e.g. the survey by Doshi( ifTTl ). Our model resembles 
the one of server with repeated vacations: a server goes 
on vacation again and again until it finds the queue non- 
empty. To the best of our knowledge, however, all existing 
models assume that the vacations are identically distributed 
whereas our setting applies to inhomogeneous vacations and 
can accommodate the case when the duration of a vacation 
increases in the average if the queue is found empty. 

We exploit the well known tool "Stochastic decomposition" 
( El) to derive the main results. Stochastic decomposition 
property of M/G/l type queueing system with server vacation 
is one of the most remarkable results shown by 1121 . The 
stationary queue length distribution at a random point in 
time can be decomposed in two or more parts where one 
part corresponds to stationary queue length distribution of 
M/G/l system without vacation. This type of decomposition 
was first observed by iTPTl . and subsequently by lfl4l . 1151 . 
EGD, H3, DID, ED- Considerable attention had been paid 
to the steady state analysis under appropriate conditions on 
vacation sequences. Most of the references can be found in 
two excellent review articles by [fTTl and fl20l in (1986). 

Our model differs from the vacation model of lfl2l due to 
the presence of inhomogeneous repeated vacation and warm- 
up time and vacation trigger time. However, the decomposition 
property is still applicable to our model since it holds the 
required assumptions stated as in ED : 

• The sequence of service times being independent of 
arrival process and independent of sequence of vacation 
periods that precede that service time. The service times 
forms i.i.d. sequence. 

> Exhaustive service and the system is in queue is stable, 
i.e. the server utilization is less than one. More over 
packets do not balk, defect or renege from the system. 

> Packets are served in First Come First Serve (FCFS) (or 
any order) which is not dependent on their service times. 



Symbol 


Meaning 


N 


Initial Queue size at the beginning of busy period , 


X(t) 


Queue size at time t , 


I 


Idle duration in one cycle , 


B 


Busy duration in one cycle , 


Vi 


Size of ith Vacation , 


c 


Expected number of vacations , 


Qn 


Amount of work in queue with N packets , 




Laplace of fcth vacation with parameter s , 


1 w 


Warm up duration (before busy period), 




Conditional warm up duration = T w l{tf > Tt}, 


N(t) 


Request arrival during time t , 


Tt 


Vacation trigger time , 


tf 


Time of first arrival if there is an arrival during Tt , 



TABLE I 

Glossary : Main notation used throughout the paper 



• Service times is non preemptive. 

• Lack of Anticipation Assumption (LAA) for vacation 
termination. 

Stochastic decomposition property allows us to obtain various 
distribution using Probability generating function (pgf) and 
Laplace in section [TTTb through highly simplified approach 
which in turns yields more insight of the system. The rest 
of the paper is organized as follows. Section [TT] describes 
our system model whose analysis is presented in Sect. [HI] 
Our modeling framework is applied to the power saving 
mechanism in a WiMAX standard through four scenarios 
in Sect. [TV] Section [V] formulates several performance and 
optimization problems whose results are shown and discussed 
in Sect.[vTJ Section [vTll concludes the paper and outlines some 
perspectives. 

II. System Model and Notation 

Consider an M/G/l queue in which the server goes on 
vacation for a predefined period once the queue is observed 
empty for a vacation trigger duration. At the end of a vacation 
period, a new vacation initiates as long as no request awaits 
in the queue. We consider the exhaustive service regime, i.e., 
once the server has started serving customers, it continues to 
serve the queue until the queue empties. Request arrivals are 
assumed to form a Poisson process, denoted N(t),t > 0, with 
rate A. Let a denote a generic random variable having the 
same (general) distribution as the queue service times. 

Note that the queue size at the beginning of a busy period 
impacts the duration of this busy period and is itself impacted 
by the duration of the last vacation period. Because arrivals 
are Poisson (a non-negative Levy input process would have 
been enough), the queue regenerates each time it empties and 
the cycles are i.i.d. Each regeneration cycle consists of: 

1) Vacation Trigger time; Failing any arrival during trigger 
time, denoted by T t , activate the vacation mode. How- 
ever vacation is deferred if there is an arrival during 
Tt which mimics the standard M/G/l queue. Time of 
first arrival, denoted by tf, is the idle duration , if the 
vacation is not triggered; 

2) an idle period; let / denote a generic random variable 
having the same distribution as the queue idle periods, 
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Fig. 1. Sample trajectory of the queue size during a regeneration cycle. 

a generic idle period / consists of ( vacation periods 
denoted V\ , . . . , V$ and vacation trigger time T t ; 

3) a warm-up period; it is a fixed duration denoted T w 
during which the server is warming up to start serving 
requests. T^, denotes the conditional warm up dura- 
tion when the vacation mode is triggered, i.e., T^, = 
T w \{t f > T t }; 

4) a busy period; let B denote a generic random variable 
having the same distribution as the queue busy periods. 

The distribution of Vi may depend on i, so the repeated 
vacations are not identically distributed. They are however 
assumed to be independent. 

Let X(t) denote the queue size at time t. It will be useful 
to define the following instants relatively to the beginning of 
a generic cycle (in other words, t = at the beginning of the 
generic cycle): 

• Vi refers to the end of the ith vacation period, for i ~ 
1, . . . , C; observe that the idle period ends at Vq\ we have 

Vi = E 5=i ^ and I = V C = Eti V; 

• Tjv refers to the beginning of the busy period B; we 
define N :— X(Tn) as the queue size at the beginning 
of a busy period; 

• Tj refers to the first time the queue size decreases to the 
value i (i.e. X(Ti) = i) for i = N — 1, . . . , 0; observe 
that the cycle ends at T . 

The times {Ti}i=j\r jv-i ... o delimit N sub periods in B, as 
can be seen in Fig. [TJ We can write B = Ei=i &i where 
Bi = Tj_x — Tj. The random variable N is in fact the number 
of arrivals from t = until time TV, even though all of the 
arrivals occur between Vf-i and T/y. Introduce Nj as the 
number of requests that have arrived up to time V£ (i.e. during 
period /) and Nt w as the number of arrivals during the warm- 
up period T w . Hence N = Nj + Nt w ■ Note that the queue 
size at the end of idle duration is X(I) = Nj. 

A possible trajectory of X(t) during a regeneration cycle 
is depicted in Fig. Q] where we have shown the notation 
introduced so far. 



the system response time. The gain from idling the server is 
introduced in the special case when the model is applied to 
study the power save operation in wireless technologies; see 
Sect. ED 

A. The Number of Vacations 

To compute the distribution of f, the number of vacation 
periods during an idle period, we first observe that the event 
C > i is equivalent to the event of no arrivals during Vi-i = 
EL=i Vk - Note that £ > 1 reflects that the first arrival occurred 
after the trigger wait time T t . Equivalently, the event £ = 
reflects that the at least one arrival occurred during the vacation 
trigger time T t . When £ = 0, let the arrival time of the first 
customer is denoted by tf conditioned such that tf < T t . 

Let Ax t denote the event of no arrival during [0,T f ] and 
£r t (A) denotes e _AT *. Let Ak denote the event of no arrivals 
during the period of time 14, and A c k denote the comple- 
mentary event. Denoting by Cu{s) '■= E[cxp(— sVu)] and 
j-(,s) := E[cxp(— sVi-i)] the Laplace Stieltjes Transform 
(LST) of Vk and V^_i respectively, we can readily write 



p(C = o) 



P(C = i) = 



[l{A c T \]=E[E[l{A c Tt }\T t ] 



E [1 - cxp(-ATi)] = l-£ Tt (A), 
P{A\) = E[l{^}] = E[E[l{^}|Vi] 
E [exp(-ATt)] E [1 - exp(-AVL)] 
= £ Tt (A)(l-A(A)), 



(1) 



(2) 



i-i 



and for i > 1, we have 

P(C = %>i) = P(A Tt )l[P(A k )P(A°) 

k=l 

P(C>i)\(i>i) = P(A Tt )T[P(A k ) 

k=l 

i-1 

= C Tt (X)l[L k (X) = C Tt (X)C^(X), 



(3) 



(4) 



k=l 



where we have used the fact that arrivals are Poisson with 
rate A. The product \fl= a £k{X) is defined as equal to 1 for 
any b < a. Let CT t (s) ■= exp(~sT t ). Using |@), the expected 
number of vacations in an idle period is given by 



![C] = £ ^P{C = i)=Y J P(t>i) = C Tt (A) fe(A). 

i-1 i=l 

(5) 



i=0 



III. Analysis 

This section is devoted to the analysis of the queueing sys- 
tem presented in Sect. [TT] We will characterize the distributions 
of £ and N, derive the expectations of /, N, B and X(t) 
and the second moments of / and N, and lastly compute 



B. The Idle Period 

The system goes on vacation only if the inactivity duration 
is more than the trigger time T t . Therefore, if the vacation is 
triggered the idle period is the sum of all the vacation durations 
including vacation trigger time. Otherwise, Idle period is only 
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the duration of first arrival ( idle period of standard M/G/l 
queue). The idle period is thus given by 



and, 



I = 



C 



P(N I =j,( = 0) = 



T t + ^^Jl c {C = 0} + T / l{C = 0} 

i=l 

C 

^inp^J + ^lK^O}. 



l-£ Tt (A), if i = 1 



0. 



otherwise. 



(6) 



Using the equality Y^i=i ^ = Y^Li ^jl{C _: expected 
idle period can be obtained as follows 

E[I] = E[mm(T t ,t f )]+E (f^ Vtl{( > i}\ 1{( ? 0} 
= E[(T / )]P(C = 0)+E[T t ]P(C^0) 

OO 

+ ^E[^e- A ^-)]/: Tt (A) 
= -£ c Tt (X) + T t C Tt (X) + C Tt (X)J2nVi\C^(X) (7) 

i— 1 

Note that when T t = 0, the idle period reduces to 

¥,[1] = Y^Li ^Wi] 111=1 ^fc(A), which is in congruence with 
HHl(eq.4). 

C. The Initial Queue Size Distribution in Busy Period 

The number of requests/packets waiting in the queue at the 
beginning of busy period is N = Nj + Nt^, where 2$ = 
T w (l{(tf > T t )}). The indicator function simply indicates 
that the warm up period is accounted only when sleep mode 
is triggered. Since the number of arrival during the idle period 
is independent of the arrival during the warm up period T w , 
the pgf of N is the product of pgf of Ni and Nt w ■ Therefore, 
we have the pgf of N given as 

N(z)=Ni(z) N Ta {z). (8) 

The queue size pgf during the idle period Nj(z) is given by 



The first term refers to the scenario when vacation is triggered, 
while the second term depicts when vacation is not triggered. 
In that case there is an arrival before the vacation trigger time 
T t . Therefore, in that case it is nothing but a standard M/G/l 
queue (without vacation). 

Denoting (A)) = (1 — £,T t (A)), the z transform of initial 
queue size Z(.) is given using eq. (O by 

oo 

Ni(z) = z m H N i = m) = zP(Nj = 1) 

m=0 
oo 

+ J2 z m F(Ni = m) 

m— 2 

oo 

= £^(A)z + ]TA(A(l-z))^(A)£ Tt (A).(10) 

i=l 

Since the arrival is a Poisson process, the pgf of arrival during 
the fixed warm up period T w is given as 



> T t )z*F(N T „ = i) 



i=0 



Nj(z) = J2 z% F ( Ni = i ) = Yl zi F ( Ni = *)• 



(9) 



i=0 



i=l 



The last equality above is due to the fact that at least one 
arrival is sure during idle period, i.e., P(Nj = 0) = 0. The 
number of arrivals during the entire idle period is the sum 
of arrivals during each vacation periods ViS, if there are any. 
To derive the distribution of Nj, we first compute the joint 
distribution of Ni and £, the number of vacations in an idle 
period. Observe that Nj takes value in N*. Noting that the 
vacation trigger event £ = is also possible, we can write 

p(JV> = j,C = »,C>i) 

= p(j arrivals in V u 1{A U . . . , A^}, l c {A Tt } 



+P(i / < T t )zJP{N TA = i)] 
= £ Tt (A)£r„(A(l-z))+£^ t (A). (11) 

Where the Laplace transform of the arrivals during the warm 
up period T w is N Tvj (z) = e-^ii-z) ._ C Tw (X(l - z)). 

The above equations combines to yield N(z) = 
Nj(z)Nt^(z), given by 

N(z) = (z£ c Tt (X)+^Ci(X(l-z))/^(X)L Tt (X)j 

^ Tt (X)£ T JX(l-z)) + C c Tt (X)). (12) 

Setting T t = corresponds the model presented in 11221 which 
is forced vacation scenario, while T t = oo corresponds to 
simple the M/G/l queue. 

Noting that, z transform is one of well known tool to 
obtain moments by using the relation N^(l) = E[N(N - 
1) . . . (N — i + 1)], which simply means the evaluation of ?ith 
derivative of N, denoted as N^ n \,), at z = 1, we derive 
the first, second and third derivatives of N(z) (which will 
be required in latter sections). The first derivative of N(z) is 
given by 



Ni{z)N$Hz) + N\ 1] {z)N Tiu {z). (13) 



where 



E 



exp(-AK t 



lC k (X)C Tt (X). 



N?\ Z ) = £^(A)+]T 



k=l 



£ Tt (X)XT w £ T JX(l-z)), 
" dd(X(l-z)) 

i=l 



dz 



(14) 

£^(^(^.(15) 
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From the definition of z transform we know that Nt iTi (1) = 
Ni(l) = 1 and thus N^(l) = XT w £ Tt (\). Further substi- 
tuting eq. (O in A^j 1 ^.), we have 



E[A/] 



A 



A 



^(A^E^^A) 



A[E[J]-T t £ n (A)] = AE[/] 



(16) 



where E[J] := E[7] - T t £ Tt (\). Combining all together we 
obtain the first moment, 

E[N] = N^(l) = N I (l)N^}(l) + N^\l)N T -(l) 

= XT w C Tt (\) + AE[7]. (17) 

Proceeding further for the second moment, the derivative of 
eq. ( fT3] l is 

a< 2) (z) = iv>(«)jyrg(2!) + JVT«(«)Jvi 2) («) 

+2iV^(z)ivj 1) (z). (18) 
where, using eq. (T5[ and eq. (fl4] i. we have 



Aj 2) (z) 



E 



-£jn<A)Ar t (A), 



(AT u ,) 2 J C Tt (A)£ Tl „(A(l-2)). 



(19) 



Evaluating at z = 1, we have 

oo 

4 2 } (1) = 5>[(A^) 2 fe(A)£T t (A) = A 2 E[/ a ], (20) 



^(l) = (Ar ro ) 2 £ Tt (A). 



(21) 



where E[I a ] := Y,Zi E K 2 ]^^(^)^T t (A). Thus, from eq. 
(IT8b we can write 

A< 2 >(1) = A 2 (T 2 £ Tt (A) + 2T lu £ Tt (A)E[7] + E[J„]) . (22) 

Therefore the second moment is 

E[A 2 ] = E[A] + N (2 \l) = A(£ Tt (A)T„, + E[7]) 

+A 2 (£ Tt (A)T 2 + 2T w C Tt (A)E[7] + E[J„]). (23) 

Note that expected initial queue size E[N] and its second 
moment E[A 2 ] obtained in l22l feq. 11,12] is a special case of 
eq. (fTTT i. ( 1231 when T t = (forced vacation scenario). Note 
that, setting the parameter T t — essentially means that the 
vacation must be triggered whenever the system is idle, which 
is a particular case. 

D. Queue Size Distribution 

Queue size is given by the number of requests/packets in 
the queue seen by any random arriving packet. From "Poisson 
Arrival See Time Average" (PASTA) ( (23]) the stationary 
queue size is equivalent to the number of requests waiting in 
the queue to be served left behind by a random departure. 
Further more, from l24l . it is equivalent to the queue size 



using the workload decomposition or stochastic decomposition 
approach (refer to IT2l ). The pgf of the stationary distribution 
of the number of request/packets left behind that a random 
departing customer is given by 



X(z) 



1 - N(z) 
E[N](l - z) 



X M/G/l( z )- 



(24) 



where Z(.) denotes the pgf of queue length at the beginning 
of busy period, -X"m/g/i(-) denotes the pgf of the number of 
customer left behind in a standard A//G/1 queue. In stationary 
regime, the distribution of queue length can also be given by 
X{z). From |25)(pp. 210), the pgf of a standard M/G/l queue 
is given by 



X M/G/l( z ) 



(l-p)(l-z)o-(A-Az) 
cr(A — \z) — z 



(25) 



where a(.) is the service time distribution. 

1) Expected Queue length: The moments of queue length 
can be directly obtained from X(z) by calculating its deriva- 
tives at z = 1. We double derivate eq. d24T i with respect to z 
and do some simple calculus to obtain the expectation of X(.) 
which is given by 

E[N] ((1 - z)X {2) (z) - 2X (1) (z)) 

= (1 - N(z))X^ /G/1 {z) - 2N^(z)X^ /G/1 (z) 



-N^(z)X M/G/1 (z) 



(26) 



Using the relations [1 — Z(l)] = (which is direct from eq. 
d24l )) at z = 1, we have 

E[X) = = —±1+E[X M/G/1 ] . (27) 

Substituting N^{1) from (|22) and X^ /G/1 {1) = 

E[X M/G/1 ] = p + g'fK 1 from iH (Pollaczek-Khinchin 

mean value formula(sec. 5.6)) in eq. d26l i, we finally obtain 
the expected queue length 



E[X] = 



A(T 2 £ Tt (A) + 2T w £ Tt (\)E{I] +E[I a }) 
2(T w £ Tt (X)+E[I}) 
X 2 E[a 2 } 



2(1 -P) 



(28) 



E. Busy period 

The length of the busy period, denoted by B, depends on 
the number of customers/packets are waiting at the end of the 
vacation interval. If there are N packet requests are waiting, 
the subsequent busy period will consists of N independent 
busy periods, each of which is denoted by B\ which mimics 
the single request service time as in M/G/l queue. Therefore, 
we have 

oo 

B*(s)=Y,Ni[B* 1 ( S )Y=N[B* 1 ( S )} 
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Thus the expected busy period can be given by 



Thus, we can express 



[B\ = E[JV]E[Bi 



ELY]E[ct] 
1-P 



(29) 



F. Sojourn time 

Assume the waiting time of a customer is independent of 
the part of the arrival process that occurs after the customer's 
arrival epoch, which is easy to show. Our policy, which is 
FCFS discipline, falls in this category. The waiting time of 
an arbitrary customer in a queue is exactly the number of 
customer ahead of the tagged customer in the queue under 
FCFS scheme. The number of customers left behind by the 
tagged customers is precisely the number of arrivals during 
the sojourn time (waiting + service) of the tagged customers, 
denoted its pgf by N(z). Since Poisson arrival see time 
average (PASTA, see Wolf), the pgf of number customer ahead 
of a random customer has the same pgf as N(z). Therefore we 
can express the LST of the waiting time W*(s) of a random 
customer in the queue as (from lfl2l ) 

W{S) = sE[N] W m/g/i(s)- (30) 

Where, W / ^ G y 1 (s) is the LST of waiting time of an arbitrary 
request in the queue (excluding its service time) of a standard 
M/G/l queue. From (261(1.45), we have W* M/G/l {s) = 

s- S \+\a*(s) ' w here LST of service time is given by u*(s) = 
E[e~ scr ]. Thus, we have 



W*(s) = 



A[l - N(l - s/X)} s{l-p) 



sE[N] 
k [1-N(1-s/X)} 
s-X + Xa*(s) 



s - A + Xa*(s) 



(31) 



where K = The moments of the W(.) can be obtained 

from its LST by simply evaluating its derivatives at s = 0, i.e., 
E[W* n ] = (-l) n W< n \0) as follows, (Refer appendix U for 
detailed computation), 



E[W] 



N^jl) AE[q- 2 ] 
2XE[N] + 2(1 - p) 



(32) 



G. Message Response Time 

The message response time T is defined as the time interval 
from the arrival time of an arbitrary message to the time when 
it leaves the system after the service completion. The mean 
message response time is said to be the single most important 
performance measure [24] (page 162) for the system without 
blocking. 

The response time of a message consists of the waiting time 
W and the service time a. Since the waiting time and service 
times are independent, we can express the LST of response 
time T and mean waiting time directly as follows 



T*(s) = W*(s)a*(s) =>■ E[T] = E[W] + E[cr]. 



(33) 



EfTl 



EH 



N™(1) AEcr 2 
2XE[N] + 2(1 - p) 

TlC Tt (A) + 2T w C Tt (A)E[7] + E[/ ]) 



2(r tB r Tt (A) + E[j] 



AEfi 



2(1 -p) 



Efcrl 



(34) 



The last term in bracket above is the mean response time 
of standard M/G/l queue, denoted by E[Tm/q/i], (refer 
11261 ). while the first part is the additional contribution due to 
vacation( which includes warm up and trigger time). Therefore, 
the expected sojourn time can be rewritten as 



E[T] = 



TlC Tt (A) + 2T w C Tt (A)E[J] + E[/ a ] 
2(T w £ Tt (X) +E[I}) 



SPm/g/i] (35) 



Remark 1: One can also obtain the expected time a cus- 
tomer spends in the queue using Little's formula as follows, 



E[T] 



ELY] 
A 



(36) 



On substitution of ELY] from ( [28] ), we obtain the same 
expression of (1351 . 

As the rate A — > 1/E [a] (recall that the stability condition 
enforces that AE [a] < 1), we must have P(C = 1) — > 1 (thus 
£i(A) — > 0) whatever the distribution of the vacations. There 
will then be only one vacation period in most idle periods. 
Therefore, at large input rates, the largest contribution to the 
sojourn time is expected to come from the waiting time when 
the server is active (queueing delays). 

H. Excess Waiting time 

Excess waiting time could be a performance metric of 
interest in context of QoS of delay sensitive class of traffic. In 
several applications (e.g. VoIP), it is important to observe the 
excess waiting time. A packet arriving late then certain delay 
is not useful and may be required to discard. Motivated from 
such necessity, we briefly discuss the excess waiting time. 
The main idea is to state that we can estimate the bound 
on excess waiting time by exploiting Markov inequality ( 
ll27l [sec. 5.5.1]). We obtain the bound on the probability of 
waiting time as follows using first and second moments as 
follows: 



Let us denote Mi = 5EL and Mo = E ^ W 2 - ■ Thus we have, 



Mo = 



E[W 2 



T3£ Tt (A) + E[/ c ] + 3(T 2 £ Tt (A)E[7] + Cr t (X)E[I a }) 



AEfcr 2 



W 2 (l — p) 



E[W] 



3w 2 (T w £. Tt (X)+E[I}) 
AELt 3 1 



3iu 2 (1 -/?)' 



(37) 



7 



Mo 



T*£T t (A) + E[/ c ] + i{TlC Tt (A)E[/] + Cr t WWa] 



AE[cr 2 ] 
w(l - p) 



Mi 



3u; 2 (T H ,/: Tt (A)+E[/]) 
AE[(T 3 ] 



3w 2 (l -p)' 



(38) 



Clearly, we can state that if ^ E[ ^ > 1, then Mi < M 2 ( 
however vice versa may not be true). This suggest that Mi is 
tighter than A/2 under condition ^ E ^j, > 1. The bound can 
be utilized to approximate the distribution of excess waiting 
time. For example, we can easily estimate the maximum packet 
dropping probability for limited buffer capacity scenario. We 
do not aim to derive a tight bound for excess waiting time 
distribution, rather we just illustrate that many of the derived 
analytic results can be easily extended to approximate the 
distributions of various performance entities, e.g. service delay, 
dropping probabilities. 



IV. Application to Power Saving 

The model analyzed in Sect. [Ill] can be used to study energy 
saving schemes used in wireless technologies. Consider the 
system composed of the base station, the wireless channel 
and the mobile node. When the energy saving mechanism is 
disabled, the system can be seen as an M/G/l queue; and 
when it is enabled, the system can be modeled as an M/G/l 
queue with vacations. The server goes on vacations repeatedly 
until the queue is found non-empty. This models the fact that 
the mobile node goes to sleep by turning off the radio as long 
as there are no packets destined to it. 

In practice, the mobile needs to turn on the radio to check 
for packets. The amount of time needed is called the listen 
window and is denoted 7). During a listen window, the mobile 
can be informed of any packet that has arrived before the 
listen window. Any arrival during a listen window can only be 
notified in the following listen window. To comply with this 
requirement, we will make all but the first vacation periods 
start with a listen window TJ. The last listen window is 
included in the warm-up period T w ( T w = TJ is considered 
for numerical evaluation). 

Let Si be a generic random variable representing the time 
for which a node is sleeping during the ith vacation period. 
We then have V 1 = S x and V t = T t + S t for i = 2, . . . , £. 
In this paper, we are assuming TJ to be a constant. As for 
the {Si}i, four cases will be considered as detailed further 
on. Figure |2(a)| (resp. |2(b)| i maps the state of an M/G/l 
queue (resp. an M/G/l queue with repeated vacations) to 
the possible states of a mobile node. 

A. The Energy Gain under Power Saving 

The performance metric defined in this section complements 
the ones derived in Sect. Hill but is specific to applications 
in wireless networks, and more precisely, to energy saving 
mechanisms. In this section, we will derive the gain in energy 
at a node should the power save mechanism be activated. 

Having noted the possible node states, we can distinguish 
between four possible levels of energy consumption, that are, 
from highest to lowest, 
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server state in 
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mobile node 
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(a) Normal mode of a mobile node. 
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(b) Repeated vacations to the possible states of a mobile node. 



Fig. 2. Mapping the M/G/l queue 



• Ghigh: experienced during exchanges of packets which 
includes the busy period(£>), 

• Gi^ten : experienced when checking for downlink packets 
which includes the listening periods TJ, 

• Gi ow : the lowest level observed when the mobile node is 
inactive but not in sleep state which includes the duration 

• G s iee P : the lowest level observed when the mobile node is 
in sleep state which includes the total vacation duration 

given by J^Li Vi- 
When the power save mechanism is not activated, the energy 
consumption per unit of time is Gi ow in idle periods (whose 
expectation is 1/A) and is equal to Ghigh during the busy 
periods (whose expectation is E[£?i]). The energy consumption 
rate can be written 



E, 



no sleep 



:= P Ghigh + (1 — p)Giow 



(39) 



where p = AE[er] = E[SJ/(1/A + E[Bi]) (loss free system). 

Consider now the case when the power save mechanism is 
activated. During busy periods, 



1 



leep 



E[G] 



(E[/l{f/ > T f }]G s leep - E[/l{t, < T t } 



CU] + T t E[l{t f > T t }](G ow - G sleep ) 
+E[T ; (C - l)l{t f > T t }](G listen - G sleep ) 
+E[T w l{t f > T t }]Gi isten + E[B]G h igh) (40) 

where the average cycle duration can be computed as E[C] := 
E[7 + T w l{t f > T t } +B] = E[J] + T w £ Tt (A) + E[B]. Evaluating 
the above terms we obtain 



+ 



(E[7] - E[J])C 5 leep + E[/]Clow + T m £ Tt CWen + E[5]C h i E h 

E[C] 

E[£ — l]Ti£.T t (A) (Cunen — CWp) . 

E[C] ( ' 
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where E[7] is defined in eq. Q and denote E[J] := j£r t (A). Note 
the following, 



E[Il{t f < T t }] = E[mm(Tt,tf)l{t f < T t }] 



+E 



J2 v MC>i})HC^o}]i{tf <T t } 

:=E[I] 



and 



E[Il{t f > T t }] = E[mm(Tt,t f )l{tf > T t }} 



+E 



£> 4 l{C>i} 1{C ^ 0}]l{t/ > T t } 



E[^.]£ Tt (A) 
(42) 

TtC-T t (A) 



£ Tt (A) :=E[7]-E[/]. 



(43) 



Observe that E[B]/E[C] = p = XE[cr] because we have assumed 
an unlimited queue (no overflow losses). The economy in energy per 
unit of time should a node enable its power saving mechanism is 
-Enosiecp — Asleep- The relative economy, or the energy gain is defined 
as 



G 



E, 



no sleep -E'sleep 
no sleep 



(44) 



p/E[B] 



Cslccp J~-T t (A) Cl 

o 



;[/]-E[/]-r ; (E[c]-i)) 

+(T l (E[(]-l)+T w )£ Tt (X) 



Cwgh 

Clisten 
Chigh 



+ 



A Chigh 

(45) 



We expect the battery lifetime to increase by the same factor. In 
practice C s iecp <C Chigh so that terms in multiplication with can 
be neglected. Letting T w =Ti, the lifetime gain reduces to 



G = 



p+^-pm 



(46) 

All performance metrics found so far have been derived as functions 
of 

• network parameters: such as the load p, the input rate A, and the 
first and second moments of the service time (E[ct] and E[a 2 ]); 

• physical parameters: such as the consumption rates Ci ow , Chigh 
and Cii sK n, neglecting C s i eep ; 

• combined physical and network parameters: such as the listen 
window Ti and warm-up period T w \ 

• the LSTs of the vacation periods and their first and second 
moments. 

In the following we will specify the distribution of the sleep windows 
{Si}i so as to compute explicitly {Lj(s)}i, {E[Vi]}i and {E[V; 2 ]}i. 

B. Sleep Windows are Deterministic 

We will first consider that the sleep windows {Si}i are determin- 
istic. More precisely, let 

C. _ _min{i-l,0 T i . • _ i 9 

u% — Li J- mm, " — A ; - - - ; 

where T m in is the initial sleep window size, a is a multiplicative 
factor, and I is the final sleep window exponent or equivalently the 
number of times the sleep window could be increased. We call T m m, 



a and I the protocol parameters. The LSTs of the vacations periods 
and their first and second moments can be rewritten 



E[vr] = 



exp(— T min s), i = 1 

exp(-(a min ^- 1 ' i >rmm + Ti)s), i = 2, 3, . . . , 

-'-mm i L i 

(a min{ l -l,0 Tm . n+Ti )n ] £ = 2 , 3, . . . , 



for n = 1, 2. We will study two cases so as to model type I and 
type II saving classes as defined in the IEEE 802. 16e standard (see 
Sect. Ull. 

Scenario D-I: This scenario is inspired by type I power saving 
classes (Tj. We consider a > 1 which implies that the first I + 1 sleep 
windows are all distinct. In particular, the value a — 2 is consistent 
with IEEE 802. 16e type I power saving classes. 

Scenario D-II: In order to mimic the type II power saving classes 
of the IEEE 802. 16e, we set a — 1 in this scenario. Letting a = 1 
equates the length of all sleep windows. Observe that we could have 
alternatively let I = 0; the resulting sleep windows would then be the 
same, namely Si = T m m for any i. Recall from Sect.|l]that in type II 
classes, a node may send or receive traffic during listen windows if 
the requests handling time is short enough. Hence, our model applies 
to these classes only if we assume that no request is sufficiently small 
to be served during a listen window Ti, 



C. Sleep Windows are Exponentially Distributed 

As an alternative to deterministic sleep windows, we explore in 
this section the situation when the sleep window Si is exponentially 
distributed with parameter p,i, for i = 1, 2, . . .. Similar to what was 
done in Sect. ITv^Bl we let 



E[Si] = — 



min{i-l,l}rp , 
u J- mm 7 



1,2,. 



(47) 



The LSTs of the {Vi}; and their first and second moments are given 
below. 



i = 1, 



E[Vi] 



1 + TminS ' 

exp(-sr ; ) 

1+a min{^U} Tm . ES ' 



min{z — 1,1 } 



E[V?] = 



9T ■ 

Zj± mm 



Tmin + Tl , 
1 



2,3,. 



2« 



2min{i— l,I}m2 



± mm I 



min{i— 1,1 } 



TminTi+Tf, 



2,3, 



Like in Sect. IIV-BI we consider two cases inspired by the first two 
types of IEEE 802. 16e power saving classes. 

Scenario E-I: Similarly to what is considered in scenario D- 
I, we consider multiplicative factors that are larger than 1, in other 

words, the values {/ii}i=i (+1 are different. When a > 1, the 

sleep windows increase in average over time. For T; = we can 
find closed-form expressions for all metrics derived in Sect. [Til] 

Scenario E-II: The last case considered in this paper is when the 
sleep windows are i.i.d. exponential random variables. This can be 
achieved by letting either a — 1 or / = in J47b . Hence pi = 1 /Tmm 
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for any i, The LSTs of the {Vi}i and their first and second moments 
simplify to 



or equivalently (recall ( 144b ) 



minimize E. 



E[Vi] 



1 ~\~ Train S 

exp(-sT t ) 
, 1 + T min s 



* = 1, 



i = 2,3, 



r m i n + T;, i = 2, 3, . . . 
IT 2 ■ i = 1 

^ mm 5 6 ^5 

2r min + 2T m inTl + T; 2 , i = 2, 3, 



V. Exploiting the Analytical Results 

Our model is useful for evaluating performance measures as a 
function of various network parameters (such as the input rate), and 
allows us to identify the protocol parameters that mostly impact the 
system performance. Instances of the expected system response time 
T and the expected energy gain G are provided in Sect. IVI-AI 

Beside performance evaluation, we will use our analytical model 
to solve a large range of optimization problems. Below we propose 
some optimization problems adapted to various degrees of knowledge 
on the parameters defining the traffic statistics. 

1) Direct optimization: This approach is useful when the traf- 
fic parameters information (e.g. the arrival rate) are directly 
available, or when they can be measured or estimated. An 
optimization problem can thus be formulated to maximize the 
system performance (e.g. the energy gain); see Sect. IV-AI for 
details. 

2) Average performance: Given that we know the probability 
distribution of the traffic parameters then we may obtain the 
protocol parameters that optimize the expected system perfor- 
mance. This optimization analysis is detailed in Sect. IV-BI 

3) Worst case performance: In the case where we do not 
have knowledge of even the statistical distribution of the 
network parameters, then we can formulate the worst case 
optimization problem which aims at guaranteeing the optimal 
performance under worst choice of network parameter. Though 
this is a more robust optimization approach, it yields a quite 
pessimistic selection of protocol parameters. Even if we do 
have knowledge of the statistical distribution, we may have to 
use a worst case performance in the case that that there is a 
strict bound on the value of some performance measure. The 
worst-case analysis will be further detailed in Sect. IV-CI 

We propose a multi objective formulation of the optimization prob- 
lem, where the performance objectives are the energy consumption 
(or performance measures directly related to the energy consumption) 
and the response time. We formulate the multi objective problem as 
a constrained optimization one: the energy related criterion will be 
optimized under a constraint on the expected sojourn time. When the 
traffic parameters are not directly known, two types of constraints 
on the expected sojourn time will be considered; in the first case 
the constraint is with respect to the average performance, and in the 
second case, it is on the worst case performance. 



A. Constrained Optimization Problem 

The objective is to optimize the protocol parameters defined earlier, 
namely, the initial window T m i n , the multiplicative factor a, and the 
exponent I. We define the following generic non-linear program: 



sleep , 



subject to T < Tqos 



(48b) 



maximize G, subject to T < Tq s 



(48a) 



where G is given in J46t . -E s iee P is given in d4 1 b and T, the 
system response time, is given in l !34t . The program J48t maximizes 
the energy gain, or equivalently, minimizes the expected energy 
consumption rate, conditioned on a maximum system response time 
Tqos- The value of Tq s is application-dependent; it needs to be small 
for interactive multimedia whereas larger values are acceptable for 
web traffic. 

The decision variables in the above optimization will correspond to 
one or more protocol parameters. For a given distribution of the sleep 
windows {Si}i, the expected number of vacations E[£], the expected 
idle period E[J], and subsequently the gain G and the expected energy 
consumption rate E s i eep will depend on the protocol parameters T m in, 
a and I and on the physical parameters Gow, Chigh and CWen (assumed 
fixed). 

We propose four types of applications of the mathematical pro- 
gram l l48t . 

1) In the first, denoted Vi, the decision variable is the initial 
expected sleep window T m i n . The parameters a and I are held 
fixed. 

2) The second mathematical program, denoted Vi, has as decision 
variable the multiplicative factor a whereas T m in and I are 
given. 

3) The decision variable of the third program, denoted Vz, is the 
exponent I. The parameters T m in and a are given. 

4) In the fourth program, denoted Va, all three protocol param- 
eters are optimized. The corresponding energy gain G is the 
highest that can be achieved. 

These four mathematical programs will be solved considering (i) 
deterministic and (ii) exponentially distributed sleep windows {Si}i. 
Instances are provided in Sect. IVI-BI 

B. Expectation Analysis 

Assume that the statistical distribution of the arrival process is 
known. Then we may obtain the protocol parameters that optimize the 
expected system performance. One may want to optimize either the 
expected energy consumption in power save mode or the economy of 
energy achieved by activating the power save mode. These problems 
are not equivalent as was the case in {48} since the energy consump- 
tion in normal mode itself also depends on the arrival process. 

As already mentioned, we consider two different constraints on the 
expected sojourn time corresponding to the situations in which the 
application is sensitive either to the worst case value (hard constraint) 
or the average value (soft constraint). 

Hard Constraints: Here, the application is very sensitive to the 
delay, so we need to ensure that the constraint on the expected sojourn 
time is always satisfied no matter the value of A. The problem is to 
find the protocol parameter 9 that achieves 

min e £ A p(A)£; s ieep(A,<9), subject to T(A,<9) < T QoS V A. 

(49) 

Another problem is to find the protocol parameter 9 that achieves 

max e J2 X P( A ) G ( A > 6> ). sub J ect t0 T ( A > e ) < T Q°s V A. (50) 

The problems j49\ and d50b are not equivalent because G depends 
also on -B no sleep which itself depends on A; recall {39}. Instances 
of {50} will be provided in Sect. IVI-Cl 

Soft Constraints: In this optimization problem it is assumed that 
the application is sensitive only to the expected sojourn time rather 
than to its worst case value. The objective is to find 8 that achieves 

min 9 E A p(A)£ s k P (A,9), subject to £ A p(A)T(A, 6) < T QoS . 

(51) 
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Alternatively, one may want to find 9 that achieves 

max<9]T A p(A)G(A,0), subject to £ A p(A)T(A, 9) < T QoS . 

(52) 

Instances of ( 1521 ) will be provided in Sect. IVI-Cl 

C. Worst Case Analysis 

When the actual input rate is unknown, then a worst case analysis 
can be performed to enhance the performance under the considered 
time constraint. Let 9 represent the protocol parameter(s) over which 
we optimize. 

Hard Constraints: Assume the constraint on the expected so- 
journ time has to be satisfied for any value of A. The problem then 
is to find 9 that achieves 

mine max A £ stap (A, 9), subject to T(A, 0) < T QoS V A. (53) 

In other words, we want to find the value of 6 that improves the 
worst possible energy consumption. A different problem consists of 
finding 6 that improves the worst possible gain, namely, 

max fl min A G(A, 0), subject to T(A, 0) < T QoS V A. (54) 

Observe that the worst possible gain is the one obtained when the 
traffic input rate tends to 1/E[<r]. Thus minA G(A, 9) « 0. Therefore, 
the above problem is meaningful only for a restricted range of small 
values of A for which the worst energy gain is far above 0. Instances 
of ([54j will be provided in Sect. IVI-CI 

Soft Constraints: Here, the application is not very sensitive to 
the delay, so it is acceptable that the constraint is respected by the 
average performance. The statistical distribution of the input rate, 
denoted p(A), is assumed to be known. The problem is to find 9 that 
achieves 

mine maxA E Attv (\, 9), subject to J2\ P(-MT(A, <?) < ^QoS- 

(55) 

Again, a different objective can be desired, namely to maximize the 
worst gain. Like what was mentioned in the previous section, the 
problem is meaningful only when the rate A is small. 

max e min A G(A, 9), subject to J2 a pW t {\ e ) < T Q°& ■ 

(56) 

Instances of d56b will be provided in Sect. IVI-CI 

VI. Results and Discussion 

We have performed an extensive numerical analysis to evaluate the 
performance of the system in terms of the expected system response 
time T given in d34t and the expected energy gain G given in d46l l; 
cf. Sect. IVI-AI In addition we have solved the problems V1-V4 for 
given values of the protocol parameters held fixed; cf. Sect. IVI-BI 
Instances of the problems l |50t , j52\ . j54\ and l |56l l are also provided; 
cf. Sect. IVI-CI We first consider vacation trigger time T t = 0, i.e. 
whenever the system is idle it is bound to go for at least one vacation. 
Latter we illustrate the impact of vacation trigger time depicting one 
of the cases. 

Physical and network parameters have been selected as follows: 

Gow/Ghigh = Gisten/Chigh = 0.2, E[ct] — 1, E[er 2 ] = 2, 
Ti=T w = 1, T QoS = 50/100. 

Unless otherwise specified, the protocol parameters are set to the 
default values: T m in = 2, a = 2, Z = 9 and Tt — in scenarios D-I 
and E-I, and T m in = 2, a = 1, I = and T t — in scenarios D-II 
and E-H. 

We have varied A in the interval (0, 1), T m i n in (1, 100), a in 
(1, 10), and Tt in (0, 10). The parameter Z takes integer values in the 
interval (0, 10). 



A. Performance Evaluation 

We have evaluated numerically the expected sojourn time T and 
the expected energy gain G in all four scenarios defined in Sects. ITV-Bl 
and HV-Cl varying the input rate A and the three protocol parameters 
T m in, a, and I. Our results will be presented in the following sections. 
First, we discuss the impact of each of the three parameters on the 
performance of the system in terms of T and G: impact of T m i n 
in Sect. IVI-A1I impact of a in Sect. IVI-A2I and impact of I in 
Sect. |VI-A3l Then, we comment on each of the performance metrics: 
comments on T are in Sect. IVI-A5I and comments on G are in 
Sect. IVI-A6I 

1) Impact of the initial window size T m \ n : We will first 
investigate the impact that the initial window size T m i n has on the 
performance of the system. For reasons that will be made clear later, 
this parameter is foreseen to be the most important parameter in type I 
like power saving classes (scenarios D-I and E-I) and it is the unique 
parameter in type II like power saving classes (scenarios D-II and 

E-n). 

a) Type I like power saving classes: We set a — 2, T t — 
and I = 9 in scenarios D-I and E-I. The results are graphically 
reported in Fig. [3] Figures |3(a)| and |3(b)| respectively depict the 
expected sojourn time T against the traffic input rate A and the initial 
sleep window size T m ; n when sleep windows are deterministic and 
exponentially distributed. The energy gain under the same conditions 
is depicted in Figs. |3(c)| and |3(d)| 

The size of the initial sleep window has a large impact on T 
for any value of A. More precisely, T increases linearly with an 
increasing T m - m for any A; see Figs. |3(a)| |3(b)| As for the gain G, 
it is not impacted by T m ; n , except for a small degradation at very 
small values of T m i n , hardly visible in Figs. |3(c)] and |3(d)1 

b) Type II like power saving classes: We set a = 1, T t = 
and Z = in scenarios D-II and E-II. The results are graphically 
reported in Fig. [4] Figures |4(a)| and |4(b)| respectively depict the 
expected sojourn time T against the traffic input rate A and the initial 
sleep window size T m i n when sleep windows are deterministic and 
exponentially distributed. The energy gain under the same conditions 
is depicted in Figs. |4(c)| and |4(d)| 

About the impact of Tmin on T and G, we can make similar 
observations to those made for type I like power saving classes, to 
the only exception that here the degradation of G at very small values 
of Tmin is more visible, especially in Fig. |4(c)| 

Observe that a larger T m j n yields a larger sleep time but it also 
reduces E[£] which together explains why the impact on the energy 
gain is not significant. 

2) Impact of the multiplicative factor a: The second param- 
eter used in type I like power saving classes (scenarios D-I and E-I) 
is the multiplicative factor a. In order to assess the impact of a 
on the performance of the system, we perform a numerical analysis 
in which the initial window size is T m i n = 2, the vacation trigger 
time Tt = 0, the exponent is Z = 9 and the multiplicative factor 
a is varied from 1 to 10. We evaluate the expected sojourn time 
T and the energy gain G both for deterministic (scenario D-I) and 
exponentially distributed (scenario E-I) sleep windows. We report 
the results inFig. [5] Figures |5(a)| and |5(c)| respectively depict the 
expected sojourn time T and the energy gain G against the traffic 
input rate A and the multiplicative factor a when sleep windows 
are deterministic. The results obtained when the sleep windows are 
exponentially distributed are displayed in Figs. |5(b)| and |5(d)| 

Interestingly enough, the multiplicative factor a does not impact 
the gain G. It impacts greatly T but only at very low input rates. 
Observe that T increases exponentially with an increasing a for small 
A which is reflected in Figs. |5(a)| and |5(b)| 

3) Impact of the exponent I: The third parameter used in type 
I like power saving classes (scenarios D-I and D-II) is the exponent 
Z. In order to assess the impact of the maximum sleep window size 
on the performance of the system, we perform a numerical analysis 
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(a) sojourn time, deterministic (b) sojourn time, exponential (c) energy gain, deterministic (d) energy gain, exponential sleep 
sleep windows sleep windows sleep windows windows 

Fig. 3. Impact of T m ; n on T and G in type I like power saving classes. 





(a) sojourn time, deterministic (b) sojourn time, exponential (c) energy gain, deterministic (d) energy gain, exponential sleep 
sleep windows sleep windows sleep windows windows 

Fig. 4. Impact of T m ; n on T and G in type II like power saving classes. 






(a) sojourn time, deterministic (b) sojourn time, exponential (c) energy gain, deterministic (d) energy gain, exponential sleep 
sleep windows sleep windows sleep windows windows 

Fig. 5. Impact of a on T and G with either deterministic or exponential {Si}i. 



in which the multiplicative factor is a = 2, the initial window size 
is Tmin = 2, the vacation trigger time T t = 0, and the exponent I is 
varied from to 10. We evaluate the expected sojourn time T and the 
energy gain G both for deterministic (scenario D-I) and exponentially 
distributed (scenario E-I) sleep windows. We report the results in 
Fig. [6] Figures |6(a)| and [6(c) [ respectively depict the expected sojourn 
time T and the energy gain G against the traffic input rate A and 
the exponent I when sleep windows are deterministic. The results 
obtained when the sleep windows are exponentially distributed are 
displayed in Figs. |6(b)] and |6(d")l 

Alike the multiplicative factor, the exponent I has a large impact 
on T only for a very low traffic input rate, and has no impact on 
G whatever the rate A. Observe in Fig. |6(a)| that T becomes almost 
insensitive to I beyond / = 7 (for small A). Here the initial vacation 
window T m i n is 2. We have computed T considering larger values 
of Tmin, and have observed that T saturates faster with I when the 
initial sleep window is larger. A similar behavior is observed in the 
exponential case for higher T; cf. Fig. |6(b)l 



4) Impact of Vacation Trigger Time T t : The fourth and the 
last parameter used in type I like power saving classes is the vacation 
trigger time T t . In order to assess the impact of the vacation trigger 
time on the performance of the system, we perform a numerical 
analysis in which the multiplicative factor is a = 2, the initial window 
size is T m in = 2 the exponent I = 9 and the vacation trigger time is 
varied from to 10. We evaluate the expected sojourn time T and the 
energy gain G both for deterministic (scenario D-I) and exponentially 
distributed (scenario E-I) sleep windows. We report the results in 
Fig. |7] Figures |7(a)| and |7 (c)| respectively depict the expected sojourn 
time T and the energy gain G against the traffic input rate A and the 
vacation trigger time Tt when sleep windows are deterministic. The 
results obtained when the sleep windows are exponentially distributed 
are displayed in Figs. |7(d)| and |7(b)| 

As expected, decreasing vacation trigger time enhances the prob- 
ability of the system to go on vacation resulting larger response time 
T and larger gain G as well for any A. 

5) The expected sojourn time T: The numerical results of the 
expected sojourn time T are reported in Figs.[3}{6] parts (a) and (b). 



12 






(a) sojourn time, deterministic (b) sojourn time, exponential (c) energy gain, deterministic (d) energy gain, exponential sleep 
sleep windows sleep windows sleep windows windows 

Fig. 6. Impact of I on T and G with either deterministic or exponential {Si}i. 





(a) sojourn time, deterministic (b) sojourn time, exponential (c) energy gain, deterministic (d) energy gain, exponential sleep 
sleep windows sleep windows sleep windows windows 

Fig. 7. Impact of Tt on T and G with either deterministic or exponential {Si};. 



As already mentioned, T is fairly insensitive to parameters I and 
a except for very small values of A. However, T increases linearly 
as Tmin increases. In scenarios D-I, E-I and E-II, as A increases, 
T first decreases rapidly then becomes fairly insensitive to A up to 
a certain point beyond which T increases abruptly. This can easily 
be explained. The sojourn time is essentially composed of two main 
components: the delay incurred by the vacations of the server and the 
queueing delay once the server is active. As the input rate increases, 
the first component decreases while the second one increases. For 
moderate values of A, both components balance each other yielding 
a fairly insensitive sojourn time. The large value of T at small A is 
mainly due to the ratio E[7 ]/E[J] (recall l|34|l), whereas the abrupt 

increase in T at large A is due to the term ^f} a I , which is the 

waiting time in the M/G/l queue without vacation. 

The situation in scenario D-II is different in that T is not large at 
small input rates A. Recall that in this scenario, all sleep window are 
equal to a constant T m i n . As a consequence, the delay incurred by 
the vacations of the server is not as large as in the other scenarios. 
The balance between the two main components of the sojourn time 
stretches down to small values of A. 

6) The expected energy gain G: The numerical results of the 
expected energy gain G are reported in Figs.[3j{6] parts (c) and (d). 
As already mentioned, G is insensitive to parameters / and a for any 

A, and sensitive to T m i n up to a certain initial sleep window size. 
The expected energy gain G decreases monotonically as A in- 
creases which can be explained as follows. The larger the input traffic 
rate A, the shorter we expect the idle time to be and hence the smaller 
the gain. 

B. Constrained Optimization Problem 

We have solved the constrained optimization program introduced 
in Sect. IV-AI as follows 

• V\ for T^ in when a — 2 and 1 = 9 (default values) with 
TqoS = 50 for scenario D-I and Tq s = 100 for scenario E-I, 



and when a = 1 or I = with Tq s = 50 for scenario D-II and 
Tqos = 100 for scenario E-II; 

• Vi for a* with T m i n = 2 and / = 9 (default values) with 
TqoS = 50 for scenario D-I and Tq s = 100 for scenario E-I; 

• Vz for I* when T m i n = 2 and a = 2 (default values) with 
TqoS = 50 for scenario D-I and Tq s = 100 for scenario E-I; 

• Va for (r m i n , a, I)* with Tq s = 50 for deterministic sleep 
windows and Tq s = 100 for exponential sleep windows. 

The optimal gain achieved by the four programs Va-Va and the 
gain obtained when using the default values are illustrated in Fig. [8] 
against the input rate A, for deterministic (Figs. [8(a)! and [8(b)] l and 
exponential (Figs. [8(c)1 and [8(d)) sleep windows. The right-hand-side 
graphs depict the optimal gain (returned by program Vi when a = 1) 
and the gain achieved under the default protocol parameter (T m i n = 
2). The most relevant observation to be made on each of Figs. |8(a)| 
and |8(c)| is the match between the curve labeled "optimal gain" (result 
of program Va) and the curve labeled "gain with T^ ilL " (result of 
program Vi). The interest of this observation comes from the fact 
that Va involves a multivariate optimization whereas Vi is a much 
simpler single variate program. The explanation for this match is as 
follows. The program Vi is being solved for the optimal T m i n . It thus 
quickly reduces the number of vacations E[£] to 1 (refer to Fig. [9) 
and thereby makes the role of both a and / insignificant. Hence, the 
energy gain maximized by Vi tends to the optimal gain returned by 
Vi. 

Comparing the optimal values of T m i n as returned by programs 
Vi and Va in the deterministic case (cf. Fig. [9] When maximizing 
the gain by optimizing T m m (program V\\ see Fig. [9(b)}, we observe 
in all scenarios but scenario D-II that, optimally, T m i n should first 
increase with the input rate A then decrease with increasing A for 
large values of A. This observation is rather counter-intuitive and we 
do not have an explanation for it at the moment. Our intuition that 
T m in should decrease as A increases is confirmed only in scenario 
D-II. 
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Fig. 8. Maximized/default gain versus the input rate A. 



Looking at the expected number of vacations E[£], should the 
optimal value T^ in be used, it appears that E[f] decreases asymptot- 
ically to 1 as A increases; see Fig. [9] The reason behind this is the 
energy consumption during listen windows and warm-up periods. To 
maximize the energy gain, one could minimize the factor multiplying 
C s ieep, in other words minimize E[£]. As a consequence, if T m i n is 
optimally selected, then the initial sleep window will be set large 
enough so that the server will rarely go for a second vacation period, 
thereby eliminating the unnecessary energy consumption incurred 
by potential subsequent listen windows. As a consequence, the 
multiplicative factor a and the exponent I will have a negligible effect 
on the performance of the system. 
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Fig. 9. Optimal parameters obtained from "P\—Tz and set to obtain E[£] 
versus A. 



C. Expectation and Worst Case Analysis 

In this section, we report the results of an expectation and a worst 
case analysis, considering the expected energy gain as performance 
metric. We will solve the problems stated in d50b . ( I52t . | |54I > and 
J56l >. The decision variable is the initial sleep window size T m i n . Each 
problem is solved for each of the four scenarios defined in Sects. HV-Bl 
and IIV-CI We consider a = 2 and I = 9 in scenarios D-I and E-I. 
Recall that we necessarily have a — 1 and I — in scenarios D-II 
and E-II. We consider that A may take five different values. These 
values and the corresponding probabilities p(A) are given in Table HT1 
The values of the parameter T m i n found for each of the problems are 
reported in Table Hill As expected, the soft constraint allows larger 



TABLE II 
Distribution of the input rate A 



A 


0.02 


0.05 


0.1 


0.2 


0.5 


pW 


0.3125 


0.3125 


0.1875 


0.1250 


0.0625 



value of Tmin as compared to the case of hard constraint. Since 
the vacation periods are of same duration for type II (D-II and E- 
II) policies, the difference between soft constraint to hard constraint 
cases are squeezed. We note that the problem of "Expected analysis" 



TABLE III 

EXPECTATION/WORST-CASE ANALYSIS: VALUE OF T min (IN NUMBER OF 
FRAMES) 





Expectation analysis 


Worst-case analysis 


Scenario 


hard constraint 


soft constraint 


hard constraint 


soft constraint 


D-I 


65 


92 


64 


64 


D-II 


96 


97 


94 


94 


E-I 


22 


50 


21 


21 


E-II 


69 


79 


62 


62 



minimizes the average of sleep gain while the "worst case analysis" 
solves the Maxmin of the Sleep gain. Obviously, the later case will 
have less variability which is reflected by the almost no difference 
of T m i n for their hard constraint and soft constraint case results. 

VII. Conclusion and Perspectives 

In this paper, we have analyzed the M/G/l queue with repeated 
inhomogeneous vacations. In all prior work, repeated vacations are 
assumed to be i.i.d., whereas in our model the duration of a repeated 
vacation can come from an entirely different distribution. Using 
transform-based analysis, we have derived various performance mea- 
sures of interest such as the expected system response time and the 
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gain from idling the server. We have applied the model to study the 
problem of power saving for mobile devices. The impact of the power 
saving strategy on the network performance is easily studied using 
our analysis. We have formulated various constrained optimization 
problems aimed at determining optimal parameter settings. We have 
performed an extensive numerical analysis to illustrate our results, 
considering four different strategies of power saving having either 
deterministic or exponentially distributed sleep durations. We have 
found that the parameter that most impacts the performance is the 
initial sleep window size. Hence, optimizing this parameter solely is 
enough to achieve quasi-optimal energy gain. 

In this paper, we have focused on deriving the expected sojourn 
time. However, it is possible to derive stronger results in means of the 
distribution of the sojourn time using the decomposition properties 
obtained in 1121 and the distributional form of Little's law 1281 . 
The queue length decomposition property 1121 states that the queue 
length in an M/G/l queue with vacations at an arbitrary epoch (i.e. 
in stationary regime) is distributed as the independent sum of (i) the 
queue length in the corresponding M/G/l queue without vacation 
and (ii) the queue length in the M/G/l queue with vacations at an 
arbitrary epoch during a non-busy interval. Given that our vacations 
are inhomogeneous, a significant portion of the derivations shall need 
to be repeated. However, we think it is worthy to investigate this 
approach and plan to do so in the near future. 

Other important research directions are considered. Namely, 

a) Other traffic profiles: It is interesting to consider more 
bursty real time traffic as well as TCP traffic. We expect that much 
of this work may have to be performed through simulations as the 
queueing analysis may become intractable. It is important to examine 
how our optimized parameters perform when a new type of traffic is 
introduced, and whether our robust design for the worst case Poisson 
traffic maintains its robustness beyond the Poisson arrival processes. 

b) Extensions of the protocol: So far our analysis enabled us 
to optimize parameters of the protocol. It is of interest to go beyond 
the optimization and to examine extensions or improvements of the 
protocol that would require to extend the theoretical framework as 
well. In particular we intend to examine rendering T m i n dynamic, by 
choosing its value at the nth idle time as a function of the Vq (or of 
its expectation) in the (n — l)-th idle time. 
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Appendix 

A. Computation of Initial queue length distribution 

Denoting (1— £.T t (A)) by £y (A)), the z transform of initial queue 
size Z(.) is given using eq. l[9} by 

Nj(z) 



= zm HNi =m) = z¥(Ni = 1) + z m HN 



i = m) 



= zC c Tt {\)+ £ 2 m 5Z E f ex P(^ AVi 
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OO 

= C c Tt (X)z + J2 E[exp(-AVi) exp(AzVS)]£jrr(A)£r t (A) 

oo 

= £^(A) 2 + ^£ l (A(l-^))£ irT (A)£ Tt (A). (57) 

i=l 

Since the arrival is a Poisson process, the pgf of arrival during the 
fixed warm up period T w is given as 



Ntv (z) = zi H N T™ = *) 

i=0 

oo 

= £[Pfo > T t )/P(iV Tfii = i) 

i=0 

+p(*/ < r t )z i p(Ar T - = ;)] 

jV T -(2) = £ Tt (A)£ Tra (A(l-^)) + £T t (A) 



(58) 



(59) 
(60) 



Where the Laplace transform of the arrivals during the warm up 
period T w is given is 

Nt w (z) = e - XT ^- z) :=£ T „(A(l-z)). (61) 
Combining both the above we can express the pgf of Z 

N(z) = N I (z)N T ^(z) 



zC c Tt {\) +X>(A(1 - z))£ irT (A)£ Tt (A) 



[£ Tt (A)£ T „(A(l- Z ))+£^(A)] 



(62) 



Note that Tt — correspond to have always vacation while Tt — oo 
forces no vacation. At Tt = 0, we have £r t (A) = e~ ATt = 1. We 
thus obtain N(z)\ Tt =o = E<=i« ^ i (A(l- z ))>C irT (A)] [£t„ (A(l- 
«))]. This is in congruence with the earlier result obtained in 1291 . 
However at Tt = oo, we have £t 4 (A) = e~ XTt = 0. We thus obtain 
N(z)\T t =oo = [z], which is again true as there is only one arrival as 
in standard M/G/l queue. 



B. Computation of Sojourn time 

We assume (general assumption) that the waiting time of a cus- 
tomer is independent of the part of the arrival process that occurs after 
the customer's arrival epoch. Our policy, which is FCFS discipline, 
falls in this category. The waiting time of an arbitrary customer in a 
queue is exactly the number of customer ahead of the tagged customer 
in the queue under FCFS scheme. The number of customers left 
behind by the tagged customers is precisely the number of arrivals 
during the sojourn time (waiting + service) of the tagged customers, 
denoted its pgf by N(z). Since Poisson arrival see time average 
(PASTA, see Wolf), the pgf of number customer ahead of a random 
customer has the same pgf as N(z). Therefore we can express the 
LST of the waiting time W*(s) of a random customer in the queue 
as (from 1121 ) 



W*(s) 



A[l-JV(l-s/A)] „,. , 

^ WWW- (63) 



Where, /a/i( s ) ' s tne LST of waiting time of an arbitrary request 
in the queue (excluding its service time) of a standard M/G/l queue. 
From (251(1.45), we have W^/a/ii 8 ) = 



s(l-P) 



M/G/iy/ s _a + Act*(s) 

of service time is given by <r*(s) = E[e~ SCT ]. Thus, we have 



where LST 
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\[l-N(l-8/\)] 8(1 -p) 

sE[N] s-X + Xa*{s) 

[l-N(l-s/\)] 



K 



s — A + Act*(s) 



(64) 



where K = ^fj^fp- The moments of the W(.) can be obtained from 
its LST by simply evaluating its derivatives at s = 0, i.e., £[W* n ] = 
( — l) n W*( n \0). The Expected waiting time is the first moment, 
given by E[W] = -W* (1) (0). This can be computed from ((64} by a 
routine but tedious calculation (two applications of L'Hospital's rule 
is required), 

Vr (1) (s)(s - A + Ao-(s)) + W*(s)(l + Aa (1) (s)) 

= K/X{N w (l-s/X)}, 
W* (2) (s)(s - A + Xa(s)) + 2W* (1) 0)(1 + Aa (1) (s)) 
+W*(s)Xa {2 \s) 



K/X 2 [N {2) (1- s/A)]. 



(65) 



Evaluating the above at s = and using the following: i) W*(0) = 
<x(0) = 1 (from the definition of LST); ii) 1 + Acr (1) (0) = 1-/9; 
and, iii) K = %pp; we obtain 

2W/* (1) (0)(1 - p) + Aa< 2 >(0) = ~ { -^- (V 2) (l)) , 



[N] 

11 W 2AE[iV] + 2(1 - p) 



(66) 



substituting iV (2) (l) from eq. (|22j and E[N] eq. 05J we obtain the 
mean waiting time as 

mw] = r 2 £T t (A) + 2r m £ rt (A)E[J]+E[J a ] AE[cr 2 ] 
2(T m £ Tt (A)+E[/]) 2(1-/0)- 

For the second moment, we perform the derivative one more time, 

H/* (3) (s)(s - A + Atr(s)) + 3VT (2) (s)(l + Aa (1) (s)) 
+3Vy* (1) (s)Aa (2) (s) + W*(s)Xa {3) (s) = K/X 3 [7V (3) (1 - s/A)] .(68) 
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Evaluating at s = we obtain the second moment, 

[ J 1 ' 3A 2 E[iV] + (1 - p) + 3(1 - p) ( ' 



Substituting N^(l) from eq. 177t . finally we have 
E[VK 



T 3 £ Tt (A) + E[J C ] + 3(7^£ Tt (A)EjJ] + £ Tt (A)E[J a ]) 
3(T w £ Tt (A) + E[7]) 



C. r/i;>c/ moment of Initial Queue distribution and sojourn 
time 

Moving further, the third derivative of the N(.) can be obtained 
similarly from eq. {18) , 

iV (3) (z) = Nj (z)N<& (z) + N<i 3) (z)N T ^ (z) + 3[N^ (z)N& {z) 

+Nf\z)N^(z)}. (71) 

iV (3) (l) = Nj(l)N^(l) + N<t 3) (l)N T€ ,(l) + 3[NP(1)NW(1) 

+iVf ) (l)iV^(l)]. (72) 

Using eq. J I Hi t we can obtain 

oo 

N ?H*) = E E [( A ^) 3 ]^(A(l- 2 ))£-r(A)£T t (A).(73) 

i=l 

oo 

iVf'(l) = E E [( A ^) :i ]^(A)£T t (A) = A 3 E[/ c ]. (74) 

!=1 

Where we denote E[J C ] := E[V?]£— T (A)£ Tt (A). Using eq. <20), 
we obtain 

A^(» = (AT„) 3 £ Tt (A)£ Tra (A(l- Z )). (75) 
N^(l) = (AT„) 3 £ Tt (A) (76) 
Combining from eq. l !73t - l !76t . we can express the 

jV< 3) (l) 

= A 3 (T 3 £ Tt (A) + E[/„] + 3[T^£ Tt (A)E[J] + E[/ a ]£ Tt (A)) (.77) 

The third moment is given by the relation E[N 3 ] = A r(3) (l) + 
3E[iV 2 ] — 2E[JV], The third moment of the sojourn time can be 
obtained similarly as above 

WW 2 ] - N(3)W + XE {W]E[a 2 ] AE[a 3 ] 

ElW ] - 3^E[N] + + 3(W) (?8) 



nw 2, = T 3 £ Tt (A) + EjJc] + 3(7^£ Tt (A)E[I] + £ Tf (A)E[J a ]) 

3(7V£ Tt (A) + E[J]) 



